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Abstract 

It is shown that the system of two coupled Korteweg-de Vries equations passes the 
Painleve test for integrability in nine distinct cases of its coefficients. The integrability 
of eight cases is verified by direct construction of Lax pairs, whereas for one case it 
£C) ' remains unknown. 

> ' 

in 

^ ! 1 Introduction 

Recently, Karasu [[L| proposed a Painleve classification of coupled KdV equations. More 
recently, we found that the classification of Karasu missed at least two systems which 
possessed the Painleve property and Lax pairs p. In the present paper, we give our 
version of the singularity analysis of coupled KdV equations. 

We study the following class of nonlinear systems of partial differential equations: 

u xxx + auu x + bvu x + cuv x + dvv x + mu t + nv t = 0, 
v xxx + euu x + fvu x + guv x + hvv x + pu t + qv t = 0, 

where a, b, c, d, e, f, g, h, m, n, p and q are arbitrary constants, and 

mq ^ np. (2) 

The condition (|2|) allows every system (Q) to be resolved with respect to ut and Vt, therefore 
our class (|l]) coincides with the class of all "nondegenerate KdV systems" studied by 
Karasu 0]. Since the "degenerate KdV systems" Q are in fact some over-determined 
systems Q reducible to lower-order ones, we do not classify them as coupled Korteweg-de 
Vries equations. 

The system (||) is "soft" in the sense that its coefficients can be changed by the trans- 
formation 



v 



' = yiu + y 2 v, v = y 3 u + y±v, t' = y 5 t, (3) 
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where yi, y 2 , ys, Va and y 5 are constants, yiy 4 ^ y 2 yz, and y 5 / 0. At every stage of 
our analysis, we use the transformation (||) in order to simplify the system (||) as far as 
possible. 

We perform the singularity analysis of the class (|l|) in accordance with the Weiss- 
Kruskal algorithm, which is based on the Weiss- Tabor-Carnevale singular expansions H, 
Ward's requirement not to analyze solutions at their characteristics ||, Kruskal's sim- 
plifying representation of singularity manifolds and which follows step by step the 
Ablowitz-Ramani-Segur algorithm for ordinary differential equations Q. Since the Weiss- 
Kruskal algorithm is well known and widely used, we omit all unessential details of our 
computations. 

It should be stressed, however, that the Weiss-Kruskal algorithm checks only necessary 
conditions for an equation to possess the Painleve property. Moreover, as explained in 
Sections |2.1.1| and 2.1.4, our analysis could miss some systems with the Painleve property. 



For these reasons, our work does not claim to be a "Painleve classification", and we prefer 
to think that it is simply a search for those systems ([!]) whose analytic properties are 
similar to such of the Hirota-Satsuma system || studied in |J. 

In Section [2|, we study the analytic properties of the systems (|l|) and find, up to 
arbitrary transformations (||), the following nine coupled KdV equations of Hirota-Satsuma 
type: 

ut = u xxx — Yluu x + 24kvv x , vt = —2v xxx + 12uv x , k = —1, 0, 1, (i 
u t = Uxxx — 12uu x , vt = v xxx — 6uv x , (ii 
Ut = u xxx — 12uu x , vt = ^v xxx — &vu x — I2uv x , (iii 



u t = u xxx — 9v xxx — Y2uu x + 18vu x + 18uv x , 

vt = u xxx — 5v xxx — Y2uu x + I2vu x + 6uv x + 18vv x , 

ut = u xxx + 9v xxx — \2uu x — \8vu x — 18uv x + 108vv x , 
vt = u xxx — llv xxx — \2uu x + 12vu x + A2uv x + \8vv x , 



(iv 



(v 



ut = u xxx — 12uu x , vt = v xxx — 12uv x , (vi 

ut = u xxx — 12uu x , vt = v xxx — 6vu x — 6uv x , (vii 

u t = u xxx — 12uu x , vt = —2v xxx + 12vu x + 12uv x , (viii 

u t = u xxx — 12uu x , vt = ku xxx + v xxx — 12vu x — 12uv x , k = 0, 1. (ix 

In Section ||, we find transformations between some of the systems (|)-(|ix|), and then 
prove the integrability of (|j)-([ix]) except (0) by direct construction of Lax pairs. 

In Section ||, we compare our results with the results of Karasu Q and give some 
unsolved problems. 



Coupled KdV Equations of Hirota-Satsuma Type 



257 



2 Singularity analysis 

2.1 Generic branches 
2.1.1 Leading exponents 

Starting the Weiss-Kruskal algorithm, we determine that a hypersurface 4>(x,t) = is 
non-characteristic || for the system ([!]) if cj) x ^ 0; without loss of generality we take 
4> x = 1, which simplifies computations and excludes characteristic singularity manifolds 
from our consideration. Since (|]) is a normal system ||, its general solution must contain 
six arbitrary functions of one variable. The substitution of 

n = n o (t)(/> Q + --- + u r (t)0 r+a + ---, v = v (t)^ + ■ ■ ■ + v T {t)4> r+P + • • • (4) 

with uqVq 7^ into the system ([[]) determines branches (i.e. the sets of admissible a, (3, 
uq and vo) and positions r of resonances for those branches. We require that the system 
(|l]) admits at least one singular generic branch, where a < or f3 < 0, the number of 
resonances is six, and r > for five of them. This requirement seems to be more restrictive 
than the Painleve property itself, therefore some coupled KdV equations possessing the 
Painleve property may be missed. 



Throughout Section 2.1, we analyze only the singular generic branches; all other bran- 



ches are considered in Section 2.2. Since the number of resonances is to be six, the set 
of leading terms of (|T|) must include u xxx and v xxx , and the admissible values of a and (3 
are as follow, depending on which of the other terms of (H) are leading and which of the 
coefficients are zero: 

• a = — 2 and (3 is arbitrary; 

• a = —4 and j3 = —6; 

• and the same with a T^- /?, which can be omitted due to (||). 

We reject systems (|l]) with a = —4 and f3 = —6 because of bad positions of resonances, 
r = —1,6,8,12 and ^(11 ± iVl59), and proceed to systems ([!]) with a = —2 and any 
integer (3. 

2.1.2 Systems with (3 > -2 

If a = —2 and (3 > —2, then e = in (jl]); moreover, p = if f3 > 0. Using (|3|), we make 
a = —12 in (|l]), and then find that no = 1 and 

W - W - 2) - 2/ + (3g)v = 25 ^ t . (5) 

Positions of resonances are —1, 4, 6, n, ri and r%, which correspond to the possible 
arbitrariness of cp, U4, uq, v ri , v r2 and v r3 , where r\ , r<i and rj, are three roots of 

r 3 + 3(0 - l)r 2 + (3/3 2 - 6/3 + 2 + g)r + 25 p P (t ) tv 1 = 0. (6) 

It follows from (||) that r\ + r<i + r% = 3(1 — (3). Since r±, T2 and r% are to be distinct 
non-negative integers, it is necessary that (3 < 0. Thus, we have {n, r2, r$} = {0, 1, 2} and 
p = for p = 0, {n, r 2 , r 3 } = {0, 1, 5} or {0, 2, 4} for (3 = -1, and / = \(3{{f3-l){f3-2)+g) 
from (H). 



258 



S.Yu. Sakovich 



If = 0, then r = -1,0,1,2,4 and 6, / = 0, and g = due to (|). We make 

oo oo 

b = by (||), substitute the expansions u = ^ Ui(t)(f) l ~ 2 and u = Vi{t)(j) 1 into ([]]), 

find recursion relations for U{ and nj, and then check compatibility conditions arising at 
resonances. We get c = at r = 4 and <i = n = 0atr = 6, therefore this case of (H) is 
equivalent to two non-coupled equations. 

When P = —1 and r = —1,0,1,4,5 and 6, we have / = and g = —6. Then (|3|) 

oo 

lets us make c = —2b to simplify computations. We substitute u = Yl an d 

i=0 

oo 

v = ^ into (||), and compatibility conditions at resonances give us h = p = b = 

i=0 

n = 0, d{m + 2q) = and (m + 2q)(m — q) = 0. Eliminating free parameters by (||), we 
get D and (§. 

When /3 = — 1 and r = —1,0,2,4,4 and 6, we have / = — | and g = —3, and make 
c = —2b by (||) again. Then compatibility conditions at resonances give us d = ih(4h — b), 
q = \p{b — 4h) + \m, n = 2bq, hz = pz = and h 2 = b 2 , where z = m + §6p — §hp. Here 
we have three distinct cases, which, after simplification by (|3|), become (|m|), (0) and (0). 



2.1.3 Systems with (3 — — 2 

When ol = (3 = — 2, ito and i>o are determined by the system 

12n + ouq + (b + c)u vq + di^ = 0, 
12i>o + eu 2 + (/ + g)u v + hv% = 0. 

If both uq and are some fixed constants (no resonances with r = 0), then (||) lets us 



make > — 2, but we have already considered those cases in Section 2,1.2| . The assumption 
that both «o and vq are arbitrary and independent (two resonances with r = 0) leads to a 
contradiction with (^). Therefore one of the variables no and vq should be arbitrary (one 
resonance with r = 0), and another one should be a function of it. According to (0), this is 
possible only if d = e = 0. Then, using (|3|), we make a = — 12 and c = —6, and ([<]) gives us 
h = 0, g = —12 — /, no = 1, Vi>o(i). Positions of resonances turn out to be r = —1,0,4,6, 
ri and r2, where n and T2 are two roots of 



r z -9r + U- f + bv {t) = 0. (8) 

Positions of resonances should not depend on vo(t), therefore we take 6 = 0, and (||) gives 
us T2 = 9 — r% and / = r\ — 9r± + 14, where r\ = 1, 2, 3 or 4 because < r\ < T2- 

When n = 1, compatibility conditions at resonances lead tom = n = g = 0, which 
violates (||), and (||) is not a system of evolution equations in this case. 

When 7~i = 2, compatibility conditions at resonances give us n = p = and q = m. 
Simplifying (|) by (|), we get @. 

When ri = 3, we find at resonances that n = q = 0, which is prohibited by (0). 

When ri = 4, compatibility conditions give usn = j) = g- m = 0orn = 5 + = 0, 
and this leads through to ([vTj]) and (|viii|) . 
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2.1.4 Systems with (3 < -2 

When a = —2 and = —3, we have b = c = d = h = in (|l|), make a = — 12 by @, 
and get uq = 1, / = —30 — |g, Vvo(t). Positions of resonances are r = —1,0,4,6, r\ and 
r2, where r\ + r2 = 12 and r\V2 = 47 + g. Since r\ and r2 correspond to the possible 
arbitrariness of v ri and v r2 , we have to study the following five distinct cases: r\ = 1, 2, 3, 4 
and 5. In the cases r\ = 1,2 and 3, compatibility conditions at resonances violate (|2|). 
When r\ = 4, compatibility conditions restrict arbitrary functions in expansions ([|). But 
in the case r\ = 5, when f = g = —12, we have only n = and g = mat resonances, and 
then, simplifying (|) by (|), get (g). 

When a = —2 and = —4, we have b = c = d = h = in <^). In this case, nvt is a 
leading term, therefore we have to consider n = and n/0 separately. There are seven 
distinct cases of positions of resonances when n = and six such cases when n ^ 0. In 
every case of those thirteen cases, however, compatibility conditions at resonances either 
contradict (||) or restrict arbitrary functions in (||). 

When a = —2 and /? < —5, we have b = c : = d = h = n = in (Q). There are ^ — |/3 
distinct cases of positions of resonances if /3 is odd and 1 — |/3 such cases if /3 is even. Now, 
using the Mathematica computer system [llj, we can check that, for j3 = —5, —6, . . . , —10 
and for all possible positions of resonances, compatibility conditions at resonances either 
contradict (^) or restrict arbitrary functions in (|j). This is very suggestive that in fact 
no systems ([[]) with (3 < — 3 possess the Painleve property. But we were unable to prove 
this conjecture within the Weiss-Kruskal algorithm for all [3 < —3, and therefore our list 
(|)-(0) may be incomplete. 



2.2 Non-generic branches 

The systems (|)-(|x|) admit many branches. The singular generic branches have been 
studied in Section |2.1| . All the nonsingular branches correspond to Taylor expansions 
governed by the Cauchy-Kovalevskaya theorem because the system (|l]) is written in a 
nonsingular Kovalevskaya form ||. Most of admissible singular non-generic branches of 
(|f)-(|ixl) correspond in fact to the singular generic expansions, where one or two arbitrary 
functions at resonances are taken to be zero. Therefore we have to study only the following 
singular non- generic branches: 

• a 



p = 


-2, 


"0 


= 2, 


v 2 o = 


1/k, r = 


: -2. 


,-1,3,4,6 and 8 for (@) with 


= 


-2, 


Uq 


= 2, 


v = 


3' ' 


-2, 


-1,3,4,6 and 8 for ©; 


p = 


-2, 


Uq 


= 6, 


v = 


10 „ _ 
3 ' ' — " 


-5,- 


-1,4,6,6 and 8 for ©; 


= 


-2, 


U 


= 2, 


v = 


3' ' 


-2, 


-1,4,5,6 and 6 for (0); 


= 


-2, 


U 


= 4, 


v = 


3' ' 


-4,- 


1,3,4,6 and 10 for (0). 



Compatibility conditions at all resonances of these branches turn out to be satisfied 
identically. Consequently, the systems (|)-(|ixj) have passed the Weiss-Kruskal algorithm 
well. 
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3 Integrability 



No two of the systems (|i|)-(|dJ) can be related by transformations (|3|), but not all of (§)-(|ix|) 
are distinct with respect to more general transformations. If we take (§) with k = and 
make v x = w, then we get exactly ( viii ) for u and w. If we make v x = w in (|i]), then we 
get exactly (|vi]|) for u and w. If we make v x = w in (|vj|), then we get flix) ) with k = for u 
and w. Therefore we can restrict our further study to and (|ix|) only. 

Since the systems (@)-(0) and (0) pass the Painleve test well, they can be expected to 
be integrable. Let us try to find their Lax pairs. We consider the over-determined linear 
system 



= A*, 



B^, 



(9) 



where A and B are some matrices and ^ is a column, and require that the compatibility 
condition of (g), 



At = B x - AB + BA, 
represents the system under consideration. We assume that 
A = Pu + Qv + R 



(10) 



(11) 



and B = B(u, 



v, v x 



where P, Q and R are constant matrices. Under 



this assumption, we get from ( [i0|) an explicit expression for B in terms of P, Q, R and 
a constant matrix S, as well as a set of conditions for P, Q, R and S. Then we try to 
satisfy those conditions, increasing the dimension of the matrices. (This is essentially the 
Dodd-Fordy version [|ll]] of the Wahlquist-Estabrook metho d Jl2[ .) In this way, we obtain 
the following results for the systems (0), (|i|), (pi|), (0) and (|x|): 



A (i) 



( u + a kv 

2 

v u — a 

\ 2 



\ 



(12) 



{{-2u x , u xx - 4u 2 + 8kv 2 + Aau + 8a 2 , 4kv a 



B (i) = {{~' 2u x, U xx 

16kv, -4kv x }, {iv 
—8u — 16a, 2u x }}; 



' X , —2kv xx + Akuv}, {—8u + 16a, 2u x , 
< x , —2v xx + 4uv, —2u x , u xx — 4u 2 + 8kv 2 — Aau + 8a 2 }, {16V, —4^ 



A 



(ii) 



2a 


u 





2 


—a 





V 





—a 



(13) 



B {ii) = {{-2u x -6au + 18a 3 , u xx - Au 2 + 3au x + 9a 2 u, 0}, {-8u + 18a 2 , 2u x + 6au-9a 



3 (ii) = l\-iu x -oau-t ±ao~, u xx 
0}, — 4uf — 3av x + 9a 2 u, vu 



uv x + 3anu, —9a 3 }}; 




(14) 
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B 



(iii) 



{{-2u x , 



4u 2 + 4au + 8a 2 , 0}, {-8u + 16a, 2u x , 0}, {-8v x , 4v xx - 4uv + 8av, 



o}}; 



A 



(iv) 



— v 
-2o- 

3 
2 



(7(§U-4 V ) 



(7 



(iv) 
' 4 
I 3 
2ct(u. e 



(7 



= {{6v x — 3cr(u — 3d) 
lizz — + 6u 2 + 4wv - 
3u x ) +6a 2 (u-3v)}, 



- 18a 3 , -Av 
18v 2 )}, {9v 



xx + Quv — 6?/ 

2 



27cr 2 , 6cj(m- 



(15) 

+ 2a(u x — 3v x ) + 6a 2 (u — 3v), 
3v) + 36cr 3 , -4v xx + 6uv - 6v 2 - 



{fa, 9v - 27a 2 



-6u, 



3<r(it-3u) - 18a d }}; 



A 



(ix) 




2 


-2k 



u + a 


v + r 







2 






u + a 




\ 



(16) 



B 



(ix) 



{{-2u x 



n. 



■4ii 2 + 4o-n + 8o- 2 , 0, 0}, {-8u + l6a, 2u x , 0, 0}, {- 



~2v x , ku xx -\- v xx 



8uv + Atu + Aav + lGar, -2u x , u xx - Au 2 + Aau + 8a 2 }, {8ku - 8v - 16ka + 16t, 2v x , 
—8u + 16cr, 2^^}}; where the index of A and B corresponds to the system, a and r are 
arbitrary parameters, and the cumbersome matrices B are written by rows. The invariant 



technique from [13] allows us to prove that the parameters a and r are essential, i.e. they 
cannot be eliminated by gauge transformations of *S>. 

The system (0) is the original Hirota-Satsuma equation || , and the Lax pair found in |11[ | 
corresponds to our result @ and ( |P2| ) with k = 1 up to a necessary transformation (g) 
and a gauge transformation of The system (^xj) represents the two (1 + l)-dimensional 
reductions of the (2 + l)-dimensional perturbed KdV equation fi"4[| , and the Lax pair @ 
and ( |l6| ) follows from the (2 + l)-dimensional Lax pair by reduction as well. 

In the case of (0), however, this approach leads us only to matrices A containing no 
essential parameters, at least for dimensions from 2x2 to 5x5; probably, the assump- 
tion ([ll]) is too restrictive. The method of truncation of singular expansions gives 
us a similar result: though the truncation procedure turns out to be compatible for the 
system (Q), the truncated singular expansions contain no essential parameters. Therefore 
it remains unknown whether the system (|v|) is integrable or not. 



4 Conclusion 

Let us make a brief comparison between our study of coupled KdV equations and the 
classification given in 

• The integrability of selected systems, equivalences between them, as well as non- 
generic branches were not studied in 

• Karasu considered only the case a = (5 = — 2 (in our notations of (||)) but did not 
use the condition uqVq ^ 0; therefore the results of llj could, in principle, contain all 



our systems with /3 > —2, i.e. (H)-(vlii), and should miss only (|ixf ) with j3 = —3. 
Our systems (|ii|), (0) and ( |ix|) are missed in the classification Q. 
Appropriate transformations (||) change Karasu's systems (13), (14), (15), (16), (17) 
and (19) (numbered as in 0) into our systems (jv|), (|vTl|) , flviiiD , @, (@) an d 
respectively, eliminating all free parameters. 
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• The system (20) from jij can be transformed by (|3|) into a system of two non-coupled 
equations and therefore is absent from our results. 

• The system (21) from [|J does not pass the Painleve test (even after a correction of 
misprints) . 

Below we give some unsolved problems. 

• Does the system (|v|) possess a Lax pair with an essential parameter? 

• Is there a system ([j]) with the Painleve property and (3 < — 3 in its singular generic 
branch? 

• Can some of the systems (§)-(0) and ( jix]) be related to each other by Miura and 
Backhand transformations? 
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